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I. INTRODUCTION 



Diffractive processes in hadron collisions are well de- 
scribed, with respect to the overall cross-sections, by 
Regge theory in terms of the exchange of a Pomeron with 
vacuum quantum numbers [l|. However, the nature of 
the Pomeron and its reaction mechanisms are not com- 
pletely known. A good channel attracting much atten- 
tion is the use of hard scattering to resolve the quark and 
gluon content in the Pomeron [2]- Such a parton struc- 
ture is natural in a modern QCD approach to the strongly 
interacting Pomeron. The systematic observations of 
diffractive deep inelastic scattering (DDIS) at HERA 
have increased the knowledge about the QCD Pomeron, 
providing us with the diffractive distributions of singlet 
quarks and gluons in Pomeron as well as the diffractive 
structure function Q. In hadronic collisions, we shall 
characterize an event as single diffractive if one of the 
colliding hadrons emits a Pomeron that scatters off the 
other hadron. A single diffractive reaction would there- 
fore correspond to, hi + hj — > hi-\- X, in which the intact 
hadron i emits a Pomeron [hi hi + W) which interacts 
with hadron j [hj -|- IP ^ X) leading to the referred reac- 
tion. Hadron i is detected in the final state with a large 
longitudinal momentum fraction. Hard diffractive events 
with a large momentum transfer are also characterized by 
the absence of hadronic energy in certain angular regions 
of the final state phase space (rapidity gaps). The events 
fulfilling the conditions of large rapidity gaps and a highly 
excited hadron remnant are named single diffractive in 
contrast to those in which both colliding hadrons remain 
intact as they each emit a Pomeron (central diffraction 
or double Pomeron exchange events). At high energies, 
there are important contributions from unitarization ef- 
fects to the single- Pomeron exchange cross section. These 
absorptive or unitarity corrections cause the suppression 
of any large rapidity gap process, except elastic scatter- 
ing. In the black disk limit the absorptive corrections 
may completely terminate those processes. This partially 
occurs in (anti)proton-proton collisions, where unitarity 



is nearly saturated at small impact parameters [j]. The 
multiple-Pomeron contributions depend, in general, on 
the particular hard process and it is called survival prob- 
ability factor. At the Tevatron energy, = 1.8 TeV, 
the suppression is of order 0.05-0.2 0, H, 0, IE]' whereas 
for LHC energy v^s = 14 TeV, the suppression appears 
to be 0.08-0.1 These corrections are, therefore, 

crucial for the reliability of the theoretical predictions for 
hard diffractive processes. 

We present below a calculation for diffractive produc- 
tion of heavy quarks in proton-proton collisions. The 
motivation is to produce updated theoretical estimations 
compatible with the scarce accelerator data on single 
diffractive charm and bottom hadroproduction [T^l 
and to obtain reliable predictions to the future measure- 
ments at the LHC. These predictions are quite impor- 
tant in the determination of the background processes for 
diffractive Higgs production and related reactions. The 
background for Higgs is a subject of intense debate in 
literature and is out of the scope of present work. Let us 
illustrate two representative examples: (a) for a low mass 
Higgs, Mh < 150 GeV, an interesting exclusive channel 
is the hh decay mode pp ^ p{H ^ hb)p [ll[; (b) the 
inclusive channel pp ^ [H ^ WW 1+1-^t)X [i3]. 
In both cases the heavy quark background contribution 
is sizable [Til, [l^, namely the DPE bottom production 
in (a) and Icpton contribution from heavy quark decays 
in (b). 

For the present purpose we rely on the Regge factor- 
ization and the corresponding corrections for multiple- 
Pomeron scatterings. Factorization for diffractive hard 
scattering is equivalent to the hard-scattering aspects of 
the Ingelman and Schlein model 0, where diffractive 
scattering is attributed to the exchange of a Pomeron, i.e. 
a colorless object with vacuum quantum numbers. The 
Pomeron is treated like a real particle. Thus, one con- 
siders that a diffractive electron-proton collision is due 
to an electron-Pomeron collision. Similarly, a diffractive 
proton-proton collision occurs due to a proton-Pomeron 
collision. Therefore, the diffractive hard cross sections 
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are obtained as a product of a hard-scattering coefficient, 
a known Pomeron-proton coupling, and parton densities 
in the Pomeron. The parton densities in the Pomeron 
have been systematically extracted from diffractive DIS 
measurements. In particular, the quark singlet and gluon 
content of the Pomeron is obtained from the diffractive 
structure function F^*-^' (xp, Q^). Recently, a new 
analysis of these diffractive parton distributions has been 
presented [H by the HI Collaboration in DESY-HERA. 

The paper is organized as follows. In the next section, 
we present the main formulae to compute the inclusive 
and diffractive cross sections (single and central diffrac- 
tion) for heavy flavors hadroproduction. We show the 
details concerning the parameterization for the diffrac- 
tive partons distribution in the Pomeron. In addition, 
we present the theoretical estimations for the rapidity 
gap survival probability factor. In the last section we 
present the numerical results, taking properly the exper- 
imental cuts, and perform predictions to future measure- 
ments in the CERN LHC experiment. The compatibility 
with data is analyzed and the comparison with other ap- 
proaches is considered. 



II. DIFFRACTIVE HADROPRODUCTION OF 
HEAVY FLAVORS 

Let us start by introducing the main expressions to 
compute the inclusive and diffractive cross sections for 
heavy flavor production in hadron colliders. The starting 
point is the inclusive total cross section for a process in 
which partons of two hadrons, hi and h2, interact to 
produce a heavy quark pair, hi + h2 QQ + X, at center 
of mass energy ^/s. At leading order (LO) heavy quarks 
are produced by gg fusion and qq annihilation while at 
next-to-leading order (NLO), qg + qg scattering is also 
included. At any order, the partonic cross section may 
be expressed in terms of dimensionless scaling functions 
fj^j''''' that depend only on the variable p |13|], 

Q k=0 

^ ± fS-"i») (4) . (1) 

where s is the partonic center of mass, mg is the heavy 
quark mass, fiR (pf) is the renormalization (factoriza- 
tion) scale, and p = s/Amg — 1. The cross section is 
calculated as an expansion in powers of as with fc = 
corresponding to the Born cross section at order 0(0^). 
The first correction, k — 1, corresponds to the NLO 
cross section at 0{al). The total hadronic cross section 
is obtained by convoluting the total partonic cross sec- 
tion with the parton distribution functions of the initial 




FIG. 1; (Color online) Total cross section for inclusive charm 
hadroproduction (solid line) and corresponding single diffrac- 
tive cross section without multiple-Pomeron correction (dot- 
dashed line). Available accelerator data are also shown (see 
text). 

hadrons, 

X a,j{s,m^Q,p],,p\) , (2) 

where the sum i,j = q,q,g is over all massless partons, 
xi and X2 are the hadron momentum fractions carried 
by the interacting partons and t = 4toq/s. The parton 
distribution functions, denoted by ff{xi,p), are evalu- 
ated at the factorization scale, assumed to be equal to 
the renormalization scale in our calculations. 

For sake of illustration, we calculate numerically the 
total inclusive cross sections for heavy flavors pair- 
production using the MRST set of partons fli^ . The low 
X region is particularly relevant for QQ production at 
the LHC as well as at Tevatron. For charm and bottom 
production, the gg process becomes dominant and infor- 
mation on the gluon distribution is of particular impor- 
tance. We compute the cross sections with the following 
mass and scale parameters: pc = '^-ruc (with rric = 1-5 
GeV) and ph = mt (with nib ~ 4.5 GeV). For the top 
quark case, we use pt — "^t (where rrit = 176 GeV). 
The present choice for the scales is based on the current 
phenomenology for heavy quark hadroproduction [T5| . 
The agreement with the total charm cross section data is 
fairly good. In Figs. [T] and [5] the numerical results (solid 
curves) are compared to accelerator data of charm and 
bottom hadroproduction [l6|, respectively. For the top 
cross section, we obtain atotiVs = 1.8 TeV) = 5 pb and 
c^totiVs = 14 TeV) — 2 nb. It should be stressed that 
sizable uncertainties are introduced by changing, for in- 
stance, quark masses and/or the renormalization scale. 
However, our purpose here is to estimate the diffractive 
ratios (J^ /crtot, which are less sensitive to a particular 
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choice. 

For the hard diiTractive processes we will consider the 
Ingelman-Schlein (IS) picture 0], where the Pomeron 
structure (quark and gluon content) is probed. In the 
case of single diffraction, a Pomeron is emitted by one of 
the colliding hadrons. That hadron is detected, at least 
in principle, in the final state and the remaining hadron 
scatters off the emitted Pomeron. A typical single diffrac- 
tive reaction is given by p + p — > p + QQ + X. In the IS 
approach, the single diffractive cross section is assumed 
to factorise into the total Pomeron-hadron cross section 
and the Pomeron flux factor 0]. The single diffractive 
event may then be written as (for equal hadrons) 
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where the Pomeron kineniatical variable is defined 



,0) 



/sy, where 



is the center-of-mass en- 



ergy in the Pomeron-hadron j system and ^/sij = ^/s 
the center-of-mass energy in the hadron i-hadron j sys- 
tem. The momentum transfer in the hadron i vertex is 
denoted hy ti. A similar factorization can also be ap- 
plied to double Pomeron exchange (DPE) process, where 
both colliding hadrons can in principle be detected in the 
final state. This diffractive process is also known as cen- 
tral diffraction (CD). Thus, a typical reaction would be 
p+p p+p + QQ + X , and DPE events thus are charac- 
terized by two quasi-elastic hadrons with rapidity gaps 
between them and the central heavy flavor products. The 
DPE cross section may then be written as, 

da^^ (hh^h + h + QQ + X) 



FIG. 2: (Color online) Total cross section for inclusive bottom 
hadroproduction (Solid line) and corresponding single diffrac- 
tive cross section without multiple-Pomeron correction (dot- 
dashed line). Available accelerator data are also shown (see 
text). 



Using the substitution given in Eq. ([¥]) we can proceed 
in write down the expressions for the single and central 
(DPE) diffractive cross sections for QQ production using 
the expression for the total cross section in Eq. 
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In order to obtain the corresponding expression for 
diffractive processes, one assumes that one of the 
hadrons, say hadron hi, emits a Pomeron whose partons 
interact with partons of the hadron /12. Thus the par- 
ton distribution xifi/hi{xi, ^'^) in Eq. ([2]) is replaced by 
the convolution between a putative distribution of par- 
tons in the Pomeron, /3/a/]p(/3, ^i^), and the "emission 
rate" of Pomerons by the hadron, fw/hixwjt)- The last 
quantity, /p /h (x^ , i) , is the Pomeron flux factor and its 
explicit formulation is described in terms of Regge the- 
ory. Therefore, we can rewrite the parton distribution 
as 



Similar expression holds for the DPE process, which 
reads as, 

o-?i^2(S'"^q) = dxi dx2 

i ^ = nn n J r j Xx 

^/ff/hi {^Ss^ fw/ii2 (^IP^) 



(1) /-^r" ^^(2) 



dxji 



dp I dtfTp/h^{xw, t) 



Xjp 



(3) 



and, now defining /{xw) = j^^dt fw/hi{xw,t), one ob- 
tains 

Xlfa/hiixi, fi^) = / rfa;ff /(xp) — /a/i>( — ,M^). (4) 



dx'4> r' dx 

Xl Xjp J X2 X]p 



X fi/TP I "TTJ'M^ j fj/w I ^'M^ j CT.u(s,m^,/x2).(6) 
\Xtp j \x-sp j 

Having presented the main expressions for single and 
central diffraction, in the following subsections we shortly 
present some details on the Pomeron structure function 
and Pomeron flux that will be employed in our phe- 
nomenological studies. In addition, an estimation for the 
multiple-Pomeron corrections is discussed. 



A. The Pomeron Structure Function 

In the estimates for the diffractive cross sections, we 
will consider the diffractive pdf's recently obtained by 
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the HI Collaboration at DESY-HERA ^. The Pomeron 
structure function has been modeled in terms of a light 
flavour singlet distribution consisting of u, c? and 

s quarks and anti-quarks with u — d=s = u = d — s, 
and a gluon distribution g{z). Here, z is the longitudi- 
nal momentum fraction of the parton entering the hard 
sub-process with respect to the diffractive exchange, such 
that z = (3 for the lowest order quark-parton model 
process, whereas < /? < z for higher order pro- 
cesses. The Pomeron carries vacuum quantum numbers, 
thus it is assumed that the Pomeron quark and anti- 
quark distributions are equal and flavour independent: 
Qjp = qjp = 2WJ^TP^ where is a Pomeron singlet quark 
distribution and Nf is the number of active flavours. The 
quark singlet and gluon distributions are parameterized 
at Qq with the general form. 



zh{z, QD^A.z^'^iX-zf^ exp 



0.01 



(1 



(7) 



where the last exponential factor ensures that the diffrac- 
tive pdf's vanish at z = 1. For the quark singlet distribu- 
tion, the data require the inclusion of all three parameters 
Aq, Bq and Cq in equation [71 By comparison, the gluon 
density is weakly constrained by the data, which is found 
to be insensitive to the Bg parameter. The gluon den- 
sity is thus parameterized at Qq using only the Ag and 
Cg parameters. With this parameterization, one has the 
value Qq = 1.75 GeV^ and it is referred to as the 'HI 
2006 DPDF Fit A'. It is verified that the fit procedure is 
not sensitive to the gluon pdf and a new adjust was done 
with Cg = 0. Thus, the gluon density is then a simple 
constant at the starting scale for evolution, which was 
chosen to be Qq = 2.5 GeV^ and it is referred to as the 
'HI 2006 DPDF Fit B'. 

Another important element in the calculation is the 
Pomeron flux factor, introduced in Eq. ([3]). We take the 
experimental analysis of the diffractive structure function 
Q, where the Xw dependence is parameterized using a 
flux factor motivated by Regge theory 



,-Bpi 



fw/p{xTP,t) = Atp ■ 2ap(t)-l ' 



(8) 



where the Pomeron trajectory is assumed to be linear, 
aw{t) = aff(O) -I- a'^t, and the parameters B^ and a'^p 
and their uncertainties are obtained from fits to HI FPS 
data [l7j . The normalization parameter A^p is chosen 
such that • "'° fw /-pi^TPi = 1 at Xjp = 3 • 10~^, 
where |imin| — n^pX^ / (1 — x^p) is the minimum kine- 
matically accessible value of \t\, nrip is the proton mass 
and |tcut| = 1-0 GeV^ is the limit of the measure- 
ment. Expression in Eq. ([5]) corresponds to the standard 
Pomeron flux from Regge phenomenology, based on the 
Donnachie-Landshoff model 18]. 



B. Corrections for multiple-Pomeron scattering 

In the following analysis we will consider the suppres- 
sion of the hard diffractive cross section by multiple- 
Pomeron scattering effects. This is taken into account 
through a gap survival probability factor. There has 
been large interest in the probability of rapidity gaps 
in high energy interactions to survive as they may be 
populated by secondary particles generated by rescatter- 
ing processes. This effect can be described in terms of 
screening or absorptive corrections p^ . This suppres- 
sion factor of a hard process accompanied by a rapidity 
gap depends not only on the probability of the initial 
state survive, but is sensitive to the spatial distribution 
of partons inside the incoming hadrons, and thus on the 
dynamics of the whole diffractive part of the scattering 
matrix. 

Let us introduce a short discussion on the usual proce- 
dure to determine the survival probability factor, estima- 
tion of its size and the uncertainties in its determination. 
The survival factor of a large rapidity gap (LRG) in a 
hadronic final state is the probability of a given LRG 
not be filled by debris, which originate from the soft re- 
scattering of the spectator partons and/or from the gluon 
radiation emitted by partons taking part in the hard in- 
teraction. Let A{s,b) be the amplitude of the particular 
diffractive process of interest, considered in the impact 
parameter, 5, space. Therefore, the probability that there 
is no extra inelastic interaction is 



<|5p>= 



^ cPh\A{s,h)\^ exp[-»(s,6)] 
J(Pb\A{s,b)\'' 



(9) 



where VL is the opacity (or optical density) of the in- 
teraction. This quantity can be computed using a sim- 
ple one-channel eikonal model or more involved multiple- 
channel model. The opacity ri(s, h) reaches a maximum 
in the centre of proton and becomes small in the periph- 
ery. Therefore, the survival factor depends on the spa- 
tial distribution of the constituents of the relevant sub- 
process. For instance, the spatial ^-distribution of sin- 
gle and double rapidity gap processes are assumed to be 
controlled by the slope B of the Pomeron-proton vertex, 
f3{t) (X exp(— _B|t|), and that there is no shrinkage coming 
from the Pomeron amplitude associated with the LRG 
in hard diffractive subprocesses. In addition, exp(— f2) is 
the probability that no inelastic soft interaction in the 
re-scattering eikonal chain results in inelasticity of the 
final state at energy s and inpact parameter b. 

As the ^-dependences of single and double-diffractive 
dissociation amplitudes are distinct, the survival fac- 
tor is different for each case. The survival factor de- 
creases with energy due to the growth of the opacity, 
whereas increases with the slope B. In order to illus- 
trate the calculation in Eq. let us consider the 
eikonal model (one channel) and assume a Gaussian 6- 
dependence for the scattering amplitude for a hard pro- 
cess, Ah{s, b) oc exp(— 6^/i?|^), with a constant hard ra- 
dius Rh- The opacity can be oversimplified in the form 
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^(«'^) = exp[— R^]. Here, A = aj>(0) - 1 and 
the soft slope is cx ln(s). In this simple case, Eq. 

([5]) can be analytically evaluated and then give the fol- 
lowing: <\S\'^>= RH7[RH,a{s)]/[a{s)]'^" . We use the 
■, a{s) = aos^ and 'f{a,x) is the 



notation Rh{s) 
incomplete Euler gamma function [20| . This rough calcu- 
lation shows the model dependence on the input values 
Rh, R{s) and a{s). We quote Ref. [13] for extensive 
review on the determination of LRG survival probability 
gap factor using the multi-channel eikonal models. 

Concerning the model dependence, the single chan- 
nel eikonal model considers only elastic rescatterings, 
whereas the multi channel one takes into account also 
inelastic diffractive intermediate re-scatterings. In gen- 
eral, current works in literature consider a more elaborate 
two or three channel eikonal model. These approaches 
describe correctly the relevant observables in pp and pp 
collisions as atot,ei, Bei, Osd and so on .5,i, Jj]. The cor- 
responding survival probabilities of single, double and 
central channels are not identical because each one has a 
different hard radius. The available experimental observ- 
ables which can be compared to the theoretical predic- 
tions of the survival probability factor are the hard LRG 
di-jets data obtained in the Tevatron and HERA d, |23| 
as well as diffractive hadroproduction of heavy bosons 
{W^ and Z^) in the Tevatron 21]. A direct informa- 
tion on the survival probability factor is obtained from 
the diffractive hard jets ratio mesured in Tevatron, where 

1|sp>^(^'l°i?Tw/ = 2-2 ± 0.8. The current theoretical 
predictions are in agreement with this measurement. 

For our purpose, we consider the theoretical estimates 
from Ref. § (labeled KMR), which considers a two- 
channel eikonal model that embodies pion-loop inser- 
tions in the pomeron trajectory, diffractive dissociation 
and rescattering effects. The survival probability is com- 
puted for single, central and double diffractive processes 
at several energies, assuming that the spatial distribu- 
tion in impact parameter space is driven by the slope B 
of the pomeron-proton vertex. We will consider the re- 
sults for single diffractive processes with 2B = 5.5 GeV~^ 
(slope of the electromagnetic proton form factor) and 
without A'^* excitation, which is relevant to a forward 
proton spectrometer (EPS) measurement. Thus, we have 



SD 
KMR- 



0.15, J0.09] and <]S']2> 



DPE 
KMR- 



0.08, J0.04] 



for = 1.8 TeV (Tevatron) J^i = 14 TeV (LHC)] 
There are similar theoretical estimates, as the GLM 
approach which also consider a multiple-channel 

eikonal approach. We verify that those results are con- 
sistent with each other. We quote Ref. [13] for a detailed 
comparison between the two approaches and further dis- 
cussion on model dependence. 

Eor sake of illustration, we have adjusted the gap sur- 
vival probability factor (for single diffraction) as a func- 
tion of center-of-mass energy using the following func- 
tional form < ]S']2 >SD (^y^) = a/{b + cln'^y/s), in the 
interval 10 GeV < -y/s < 14 TeV. This functional form is 
inspired in the fact that upon reaching the unitarity limit 
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FIG. 3: (Color online) Single diffractive cross section for 
charm and bottom hadroproduction considering multiple- 
Pomeron corrections as a function of energy. A comparison 
to experimental data is also shown (see text). 



TABLE I: Model predictions for single (SD) and central 
(DPE) diffractive heavy flavor production in Tevatron and 
LHC. Numbers between parenthesis represent the estimates 
using the single Pomeron exchange. 





Heavy Flavor 


i?SD (%) 


RUPB (%) 


1.96 TeV 


cc 


3.5 (23.4) 


0.08(1.2) 


1.96 TeV 


hh 


2.5(16.9) 


0.05 (0.6) 


1.96 TeV 


tt 


1.5 • 10-=^ (0.03) 




14 TeV 


cc 


2.8(30.7) 


0.08 (2.0) 


14 TeV 


hi 


2.1(23.5) 


0.035 (0.9) 


14 TeV 


tt 


0.4 (4.8) 


8 ■ 10"^ (0.23) 



the fraction of diffractive events is expected to vanish as 
1/ In y/s ^] . For energies lower than ^/s < 540 GeV the 
estimates from Ref. [3] are used. The following param- 
eters are obtained: a = 2.062, b = 4.937, c = 4.2 • lO'^ 
and d = 3.793. This is useful to obtain the single diffrac- 
tive cross sections as a function of energy when they are 
corrected by multiple-Pomeron scattering. 



III. RESULTS AND DISCUSSION 

In what follows, we present predictions for hard diffrac- 
tive production of heavy flavors based on the previous 
discussion. In the numerical calculations, we have used 
the new HI parameterizations for the diffractive pdf's 
y. The 'HI 2006 DPDF Fit A' is considered, whereas a 
replacement by 'HI 2006 DPDF Fit B' keeps the results 
unchanged. For the pdf's in the proton we have consid- 
ered the updated MRST parameterization [3]. In addi- 
tion, we have used the cut < 0.1. The single Pomeron 
results, given by Eq. ([5]), are presented in Figs. [1] and 
[5] for charm and bottom production (dot-dashed curves). 
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The single difFractive contribution is large, being of or- 
der 15-20 % from the total cross section. For the top 
case, which is not shown in a plot, we get a small cross 
section. Namely, a^f = 1.5 • lO'^ pb at = 1.96 TeV 
and af^ = 100 pb for ^/s = 14 TeV. The cross sections 
for each heavy flavor rise with energy and differ approx- 
imately by a factor proportional to I/ttiq. The results 
corrected by unitarity suppresion are a factor about 1/10 
lower than the single-Pomeron ones. In Fig. [3]we present 
the numerical results for the single diffractive production 
of charm (long-dashed curve) and bottom (solid curve) 
considering the suppression factor. We have multiplied 
Eq. ^ by the energy dependent factor < \S\'^ ("v/s), 
as discussed in previous section. Of course, the present 
curves have to be considered as average results, since the 
absolute cross sections are sensitive to different choices 
for the quark masses and renormalization scale. 

The results using the multiple-Pomeron correction can 
be somewhat compared with experimental estimates for 
the single diffractive cross section. However, data for 
diffractive production of heavy flavors are scarce. For 
charm, we quote the most accurate data for produc- 
tion of D*^ from the E690 experiment at Fermilab at 
\/s = 40 GeV [9]. The cross section is cr{pp p[cc] X) = 
0.61±0.12(stat)±0.11(sj/st) /ib, which is integrated over 
xf > 0.85. This data point is represented by the filled 
circle in Fig. [H For the bottom case, the estimation 
a{pp pMX) = 1.4 ± 0.7 /zb at = 1.8 TeV is taken 
from Ref . [23| . This value is obtained using theoretically 
predicted inclusive cross section 140 fih and the experi- 
mentally measured fraction i?^^ of diffractively produced 
beauty [l3|- This value is represented by the filled tri- 
angle in Fig. [31 Our curves reasonably describe these 
scarce data on diffractive production of heavy flavored 
hadrons. 

Concerning the DPE cross sections, we verify that the 
single-Pomeron predictions are still large. For instance, 
using Eq. © one has crj?/'^ = 58 ^ib and cr^_PE ^ i ^ 
at =1.96 TeV. The top production is suppressed 
due to kinematic threshold related to the cut-off in the 
Pomeron spectrum, Xjp < 0.1. This gives little room to 
observe diffractive top quark events at the Tevatron but is 
promising for the LHC, where the heavy flavour threshold 
suppression is less severe. The values reach cr°/^ = 460 
Atb, CT^_PE ^ and a]^jf^ = 5 pb for = 14 TeV. 

The results corrected by multiplc-Pomeron suppression 
factor are a factor about 1/100 lower than the single- 
Pomeron ones. 

Let us now compute the diffractive ratios. The sin- 
gle diffractive ratio is defined as i?sD = o-^^/o'tot and 
the central diffractive ratio by i?DPE ~ cr^^^ / <^tot- The 
results are summarized in Table I, where the diffractive 
ratios for heavy flavor production are presented for Teva- 
tron and LHC energies. The multiple-Pomeron correc- 
tion factor is taken from KMR model. The numbers 
between parenthesis represent the single Pomeron calcu- 
lation. Based on these results we verify that the charm 
and beauty production in single diffractive process could 




[GeV] 



FIG. 4: (Color online) The single diffractive ratio, 
RsD{\/s,'mQ), for charm and bottom hadroproduction con- 
sidering multiple-Pomeron corrections as a function of energy. 



be observable in Tevatron and LHC, with a diffractive 
ratio of order 2-3 %. The top production case is prob- 
lematic, mostly due to the kinematic threshold. The pre- 
dictions for central diffractive scattering are still not very 
promising. However, the study of these events is worth- 
while, at least for charm and bottom production. The 
predictions can be taken as a guide since a more detailed 
investigation with detector acceptance is deserved. For 
illustration, in Fig. |4]the energy dependence of the single 
diffractive ration for charm and bottom is presented in 
the interval 0.63 < ^/s < 14 TeV. The ratios decrease on 
energy due to the gap survival probability factor which 
behaves as 1/ In y/s. 

Our calculation can be compared to available litera- 
ture in diffractive heavy flavor production. For instance, 
in a previous calculation presented in Ref. j2^, the au- 
thor has considered the leading order approximation for 
the hard scattering and old parameterizations for the 
Pomeron structure function. Corrections for multiple- 
Pomeron exchange were not considered. Ours results for 
the absolute cross sections are not comparable but the 
values for the diffractive ratios for the single Pomeron 
case are similar (at least for charm and bottom). The 
present calculation is in agreement to the recent investi- 
gations within the light-cone dipole approach |23|. Our 
results for the single diffractive production of charm and 
bottom are similar to those in Ref. [1^, including the 
overall normalization and energy behavior. For the top 
case, our results seem to be smaller than the dipole cal- 
culations. 

In summary, we have presented predictions for diffrac- 
tive heavy flavor production at the Tevatron and the 
LHC. In calculations rely on the Regge factorization 
(single-Pomeron exchange) supplemented by gap survival 
probability factor (correction for multiple-Pomeron ex- 
change). For the Pomeron structure function, we take 
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the recent HI difFractive parton density functions ex- 
tracted from their measurement of F^^^\ The results 
are directly dependent on the quark singlet and gluon 
content of the Pomeron. We did not observe large dis- 
crepancy in using the different fit procedure for diffac- 
tive pdf's. We estimate the multiple interaction correc- 
tions taking the theoretical prediction a multiple-channel 
model (KMR) , where the gap factor decreases on energy. 
That is, < >~ 15 % for Tevatron energies going down 
to <|S'p>~ 9 % at LHC energy. The consideration of 
other models for the suppression factor does not intro- 
duce dramatic changes. We found that at the Tevatron 
single and central (DPE) diffractive charm and bottom 
quark production is observable with a single diffractive 
ratio Rg]? between 3.5 % (charm) and 2.5 % (bottom). 



The DPE cross section for charm and bottom produc- 
tion lies above the total inclusive cross section for the 
top quark and might also be observable. The diffrac- 
tive cross sections will be larger by at least one order 
of magnitude at the LHC and also diffractive top quark 
production could be observable. Therefore, the LHC can 
be a laboratory for diffractive scattering studies. 
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